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Abstract 

Monte Carlo(MC) simulations produce evidence that annealed copolymers incorporating two 
interconverting monomers, P and H, adsorb as homopolymers with an effective adsorption energy 
per monomer, e e //> that depends on the PH equilibrium constants in the bulk and at the surface. 
The cross-over exponent, <I>, is unmodified. The MC results on the overall PH ratio, the PH ratio 
at the surface and in the bulk as well as the number of adsorbed monomers are in quantitative 
agreement with this hypothesis and the theoretically derived e e //. The evidence suggests that the 
form of surface potential does not affect 3> but does influence the PH equilibrium. 
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I. INTRODUCTION 



The interest in polymer adsorption reflects both its practical importance[l] and the as- 
sociated theoretical challenges [2]. Within this domain the adsorption behavior of neutral 
water soluble polymers (NWSP), exemplified by polyethylene oxide (PEO) and polyvinyl 
pyrrolidone (P VP) [2I , occupies a special position. On the applied side this is because of their 
importance in the formulation of water-based colloidal solutions^, [|5, [f| of practical interest. 
From a theoretical point of view the adsorption of NWSP poses a distinctive problem: While 
NWSP such as PEO are homopolymers comprising of chemically identical monomers, they 
are modeled as "two-state polymers" whose monomers interconvert between hydrophilic (P) 
and hydrophobic (H) states 3, S Q, II, 12, 13, 14, 15]. These two-state models are in- 



voked in order to rationalize a phase diagram exhibiting a closed insolubility loop, with both 
upper and lower critical solution temperatures, that is thought to characterize NWSP. The 
models differ with respect to the precise identification of the two interconverting states. Nev- 
ertheless, within all of these models NWSP exhibit an annealed sequence of HP states with a 
temperature (T) dependent H/P ratio. The adsorption behavior of such annealed two-state 
polymers is expected to differ from that of homopolymers comprising of single monomer 
species and of quenched random copolymers with a fixed sequence and composition. Thus 
far, it was investigated using a self consistent field theory 16] of a laterally uniform adsorbed 
layer, comprising many adsorbed chains, in contact with a polymer solution. In the following 
we consider, in contrast, the adsorption transition of a single, terminally anchored, two-state 



polymer and compare it to the correspondin g re sults for homopolymers |2j, [17, 



and for quenched copolymers 
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25j. In particular, we investigate the adsorp- 



tion of a non-cooperative two-state polymer using Monte Carlo simulation supplemented by 
simple theory. Our discussion is mostly concerned with swollen chains under good solvent 
conditions. Simulation evidence suggests that the adsorption of both homopolymers and 
of quenched copolymers exhibits identical scaling behavior upon introduction of the appro- 
priate effective adsorption energy per copolymer's monomer, e e //. The cross-over exponent 
$ characterizing the second-order adsorption transition of these two systems is identical 



21 



24 



251 ]. With this in mind we expect similar behavior for the annealed two-state poly- 
mers. Accordingly we first identify the e e ff for the annealed copolymers case and then 
analyze the simulation data assuming that such chains adsorb as homopolymers with e e // 
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specifying the monomer-surface interaction. As we shall see, this picture is consistent with 
the simulation results. In particular, it allows to collapse the simulation data concerning 
the number of adsorbed monomers N$ onto universal curves and to reproduce the simula- 
tion results concerning the total H fraction as well as the H fraction in the bulk and at 
the surface. Importantly, in distinction to quenched copolymers, the total H fraction as 
well as the H fraction in the bulk and at the surface are not fixed. These results are of 



interest from two points of view, 
transition of homopolymersj^, 17 



irst, they complement earlier results on the adsorption 
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19 



20j and quenched copolymers [2 ll. |22. 



23 



24 



25]. 



Second, they provide a starting point for the formulation of a phenomenological free energy 
description 20J of the adsorption of NWSP modeled as two-state polymers. 

In formulating the problem we aimed to capture the generic features common to the 
various two-state models. These differ with respect to the precise identification of the inter- 
c— g Sta t^po lar v S apolarfl, heated v S nonhydratedfly QQQ, o, c,u Stere d 



vs nonclustered 



111 ] . With this in mind we focused on the simplest case where the monomers 



undergo a unimolecular HP interconversion. In the following we confront simple theory 
with off-lattice Monte Carlo simulations of the adsorption behavior of a single two-state 
polymer within this minimal model. There is no explicit solvent in the simulation and the 
monomers states are modeled as Lennard- Jones particles having identical collision diame- 
ters. The interactions between of the various monomer-monomer pairs are identical but the 
interaction parameters with the surface differ with the monomeric state. As a result the 
PH inter convesion at the surface and at the bulk are associated with different equilibrium 



constants (Figure 1). Our model is thus closest to the one proposed by Karlstrom 



of the two-state model is obtained in section II using a partition function method similar to 



that of Moghaddam and Whittington 



221 ] as well as its free energy counterpart. The conse- 



quences for the scaling analysis and blob description are discussed in section III. In section 
IV we present the details of the simulation model and the simulation results are discussed 
in section V. 
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II. THE EFFECTIVE ADSORPTION ENERGY e eff OF ANNEALED COPOLY- 
MERS 



Two methods can be used to identify e e //, the adsorption energy of an "effective" 
monomer at the surface. In one, first used by Moghaddam and Whittington to obtain bounds 



for quenched copolymers 



22 



231 ]. e e ff is obtained upon recasting the annealed copolymer 
partition function into a homopolymer-like form containing factors of the form exp(iVs^r) 
where k is the Boltzmann constant. Similarly, it is possible to consider the free energy 
of an annealed copolymer and recognize e e // in a term of the form —Nse e ff. These two 
equivalent methods yield as expected an identical e e ff. We briefly discuss the two versions 
in order to establish the relationship both to the Moghaddam- Whittington method and to 
the phenomenological free energy approach. 

We begin with the free energy of an adsorbed two-state chain. It comprises two terms, 
F con f and F seq . The first, F con f, allows for the adsorption induced confinement of the chain. 
It reflects the loss of configurational entropy of the flexible backbone upon confinement to 
a slab of thickness D ^ Rp where Rf ~ N 3 ^ 5 is the Flory radius of the swollen chain. F con f 
depends only on D, or equivalently N$, and its precise functional form is not important for 
the first part of our discussion. The second term, F seq , accounts for the standard chemical 
potentials of the P and H monomers in the bulk and at the surface as well as the mixing 
entropies at the surface and in the bulk. In the good solvent regime the chains are swollen, 
the monomer concentration is dilute and monomer-monomer interactions have negligible 
effect on the P ^ H equilibrium. The PH sequences in the bulk and at the surface are 
thus modeled as ideal one-dimensional mixtures and 

F seq (N s , x H , x SH ) = E- TS mix = N B f B + N s f s = Nf B + N s (fs - fa) ■ (1) 

Here f B and fg are respectively the free energies per monomer in the bulk and at the surface 
as specified by 

fs = xshPsh + (! - x sh)^°sp + kT [xsh^xsh + (1 - x S h) ln(l - x S h)] , ( 2 ) 
f B = xh/jPjj + (1 - x H )fi° P + kT [x H In x H + (1 - x H ) ln(l - x H )\ , (3) 

where the iV monomers in the chain comprise Ns adsorbed monomers at the surface and 
N B = N — Ns nonadsorbed monomers in the bulk. The H fraction among free monomers is 



xh = ^fjj^- while for the adsorbed monomers it is xsh = ^ff- where N BH and N SH denote 
respectively the number of H monomers in the bulk and at the surface. The standard 
chemical potentials of the different species are denoted by Note that fs and fs as 
specified by @ and ([31) are similar to the Zimm-Bragg free energy for the helix-coil transition 



in the case of zero cooperativity 



261 ] . At this point fs and fs are completely decoupled from 



each other and from F con f. The equilibrium conditions dfs/dxH = and dfs/dxsn = 
lead to the mass action laws in the bulk 



N 



BH 



X H 



N. 



K 



BP 



x H 



B 



exp 



I^h /fp 
kT 



(4) 



and at the surface 



SH 



x S h 



V , - K s = ^(- ^\/ SP )- (5) 

N SP 1 - x SH V kT J 

Accordingly, a mass action law of the form K = or x = K/(l + K) leads, upon 
substitution in (T5]) and ([3]), to the equilibrium free energies per monomer at the surface and 
in the bulk 



fs = Hsp -kT\n(l + K s ) = -kT In 



exp 



^SP 

' kT 



+ exp 



= e B = Hp-hT m(l + K B ) = -kT In 



^SH 

' kT 



»p\ ( (*H 

ex P |-- +exp -- 



(6) 



(7) 



In turn, these yield F seq (N s , x H , x S h) = Nf B q + N s (fs 9 - fs) = Ne B ~ N s e ef f at equilib- 
rium and 



-eff = Ib ~ fs =Pp- Vsp + kT\n — 



1 + K & 



+ K 



kT In 



B 



exp 



kT 



exp 



7& +exp -Sft 



• (8) 



Our sign convention for e e // follows the custom in the phenomenological theories of 
adsorption [20, bib | ■ 

Having obtained e e // from the free energy of an adsorbed chain we turn to the partition 
function argument. The partition function of an annealed two-state polymer at a surface is 

N S H=N S N BH =N B 

Z = H E E 

N S N SH =0 N BH =0 
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(N-N S )\ 



NbhKN 
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-N s 
N s \ 



N. 



BH , 



cxp 



N SH \{N S -N SH )\ 



cxp 



/4 

kT 
' kT 



Nbh 



cxp 



\ 1 N-Ns-Nbh 



exp 



-, Ns-Nsi 



(9) 



Here C^(Ns) is the number of chain trajectories with iVs surface contacts which is assumed 
to be identical to that of a homopolymer. N$p = Ns — N$h is the number of surface 
monomers in P state and Nbp = N — Ns — Nbh is the corresponding number among bulk 
monomers. Z as given by (Q counts all possible free and adsorbed PH sequences and assigns 
to each sequence the appropriate Boltzmann factor. The summations over Nbh and Nsh 
yields 



7° 

^HB 



exp 



N S 



exp 



exp 



Vsh 



— + exp — — 



kT 



+ exp 



' kT 



kT 

1 N S 



N-N s 



(10) 



^)+exp(-^ 



kT 



N 



exp ( — -£§?■) is the HP contribution to the bulk 
partition function of a two-state chain comprising N monomers. It counts the HP sequences 
in the absence of a surface or in the presence of a perfectly non-adsorbing surface. We now 
introduce C^(N S ) = C^(N S ) exp (— t§?-) which allows for both the backbone configurations 
and the HP sequence at a non-adsorbing surface. In turn, this yields a " homopolymer- like" 
Z 



N S 



cxp 



kT 



^<5+(iV s )exp(iV s %^, 



allowing to identify e e //, as given by equation (jHJ). 

The free energy formulation translates into the partition function description via 
F = — kT In Z. In particular, for each set Ns,xh,%sh the sequence partition 



seq 



function is approximately Z seq (Ns,XH,xsH) = exp — — S ^ H,XSH > an d Z_ 
En s E aH Exsh ^(^g) ex P [- F ^ {Ns kT H ' XSH) ] where C+(iV 5 ) containing the information 
on N s corresponds roughly to exp (^— Fc ^ f ) • This correspondence is incomplete because it 
is based on the Stirling approximation and the implicit assumption that 1 <C Ns < N, 
1 < N B < N, 1 < N SH < N s and 1 < N BH < N B . Alternatively, one may be- 
gin with the partition function and obtain the free energy formulation upon replacing 
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the factorials x\ in Z seq by exp(x\nx — x). As we discussed, e e // can be obtained from 
the minimized free energy which corresponds to the maximal term of Z. Minimization 
of F yields simultaneously e e ff and the equilibrium constants K p and K$- K p and Kg 
are obtainable from Z upon minimizing of the individual terms subject to the caveats 
noted earlier. The Moghaddam-Whittington(MW) annealed approximation of the adsorp- 
tion of quenched copolymers 22] provided the starting point of our discussion of the an- 
nealed copolymers adsorption. Within this approximation a quenched copolymer with a 
fixed average H fraction, N H /N = p, and a quenched sequence is modeled as a copoly- 
mer with a fixed p but with an annealed sequence 29(. The two treatments differ in 
two respects: (i) within the MW treatment /b = 0. (ii) The standard chemical po- 
tentials of monomers at the surface, fi SH and replaced in the MW approach by 
P>°sh = P-sh ~ ln P and fisp = Psp ~ ln(l -p). Accordingly F seci (Ns, x H , x S h) = N s f s where 
fs = xshPsh + (! - x sh)p° S p + kT[x S H^x S H + (1 - x SH ) ln(l - x SH )} and the equilib- 
rium condition dfs/dxsH = leads to the mass action law , X _ SH = Kg = -r^K$ and to 

-L ^ S H J- P 

equilibrium xsh = iJ^p Ks instead of the x,sh = i+k s as obtained for the fully annealed 
polymer considered by us. For the particular case considered by MW, where n SP = and 
Psn/kT = —6mw, this leads to 



f eq r / 77° \ / 77° 

z eff _ JS_i„ / PSP \ , _ / PSH 



kT kT 



exp I — — — + exp 



kT * V kT 



\ia\pexpeMW + 1 -p] , (12) 



while the partition function (Q is replaced by its MW counterpart 

N S H=N S N , 

2 = E E CUN S ) N HN S - , ] P NsH (l-p) Ns - Ns »exp(N SH e MW ). (13) 

One may thus consider the MW partition function as a special case of equation ( JTUl) in 
which n SH /kT = —euw ~ l n p, fi SP /kT = fx p /kT = — ln(l — p) and jJ^/kT = — In p. The 
two procedures also differ because the MW approach is an approximation when applied to 
quenched copolymers while its counterpart, as described above, is rigorously applicable to 
the annealed the two-state polymers considered here [30]. 



III. SCALING ANALYSIS AND ADSORPTION BLOBS 



Having identified e e // of two-state polymers we are in a position to discuss the scaling 
analysis of such polymers and the corresponding blob picture. For two-state polymers it 
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assumes, as is the case for quenched copolymers, that the polymers adsorb as homopolymers 
but with an excess adsorption energy per monomer of 



ekT 



e eff ~ 6 eff 



(14) 



Here e c e ^ is a constant, model dependent, critical adsorption energy at the limit of 



T = ( £ e// " tiff) A 



eff 



0, N —>■ oo while eiV* = const' 



17j. In contrast to simple ho- 



mopolymers where e ~ const' [20j the T dependence of two-state polymers e is 

exp 



e(T) =ln 



kT 



exp ( ) + exp ( 



ft 



kT 



C eff 
kT ' 



(15) 



Upon identifying e, the " homopolymer-like" scaling hypothesis for N$ is standard. In par- 
ticular 



N s « N*g s (x) where s = rA^*, (16) 
where <7 s (x) is a universal scaling function 17, 20]. At the transition, where Ns ~ A^*, 



we thus require g s (x = 0) ~ 1. In the adsorption region, x » 1, <? s (x) follows a power law 
behavior, x qs . Since N s ~ N when x >> 1 this leads to q s = (1 - $)/$ or A^ rs jVgC 1 -*)/*, 
Overall, a plot of Ns/N® vs rA^* should collapse the data onto a single curve corresponding 
to g s (rN^). To this end it is first necessary to determine the model dependent e c e ^. For 
$ = 1/2, the currently accepted value, 

M „ A^- 1 ~ iV°, (17) 

at the critical point, e e ff = e c e ^. This indicates that curves of N^/N vs. e e ff for different N 
values intersect at e e ff = e c e ^ thus allowing to determine e c e ^ from the intersection point [31I . 

We now return to the free energy per adsorbed chain, F. It comprises of two terms 
F = F con f + F seq . The adsorption free energy F seq ps —NsekT allows for the attractive 
monomer-surface contacts. F con f reflects the confinement of the polymer into a slab of 



thickness D < Rp. Within the blob picture 



20fl F conf « kTN b iob ~ kTf where N blob is the 
number of confinement "D blobs" comprising each of g monomers such that g v a ~ D. Since 
N s ^g® we obtain g ps (^) * _1 and thus 

F(N S ) fN s \ T ^ , . 

"ir 1 ^ bf < 18 > 



minimization with respect to Ng than yields the equilibrium Ns 

N s ~ Ne 1 ^. (19) 



and the equilibrium free energy of the adsorbed chain 

- -Net. (20) 

In turn, the equilibrium Ng together with xh = K B /{1 + K B ) and xsh = Kg/(1 + Ks) 
determine the total number of H monomers N B = Ng l +f Cs + (N — Ng) X ^§ B or 

N H i=* if s / if B /olN 

" e * * 77^- ( 21 ) 



AT 1 + V / 1 + K B 

As noted earlier, the recent consensus regarding the cross-over exponent suggests $ = 1/2 
The corresponding adsorption blobs for three dimensional, self-avoiding chains are 
specified by 

g fa 1/e 2 ; D « a/e 6/5 ; N 5 « JVe. (22) 

However, simulation results suggest that finite size effects are important and for finite chains 
$ ~ 3/5 is more realistic [3]. For relatively short chains the appropriate blob picture may 
thus involve 

g w l/e 5/3 ; D fa a/e ;N S « Ne 2/3 . (23) 

For the adsorption of free chains, the form of e e // leads to distinctive adsorption constant 
in dilute surface regime, when the adsorbed chains do not overlap. In this regime the 
chemical potential of the adsorbed chains is fi a d s ~ F eq + kT In V where V is the activity of 
the adsorbed chains and F eq fa —kTNe^ is the standard chemical potential of an adsorbed 
chain. The chemical potential of the free chains in the bulk is Hbuik ~ kT In cut where c^uik 
is the activity of the free chains. The adsorption isotherm, as obtained from fi aa > s = fibuik, 



is T fa CbuikKads where K a d s ~ exp(A/e*) is the adsorption constant 20]. K a d s for two state 
polymers and $ = 1/2 assumes thus the form 



Kads = expN < In 



<"-xp [ ~^§f) + exp 



expl-£f) +exp [~W 




(24) 
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The T dependence of K ads as given by ([21]) differs from that of K ads = exp (iV^ 1 /*) as 
obtained for homopolymer adsorption with e = 5 ~ const' [201] . 

As we shall see in section IV the approach discussed above accounts well for the simulation 
results. Confrontation with experimental results is more difficult. The four parameters 
involved, fi° H , /i P , fj,° SP and n° SH are specific both to the particular NWSP considered and 
to the model used in order to analyze the data. While fi° H , /z P were determined for a 
number of two-state models, the full set of parameters was only determined for PEO within 
the Karlstrom model by fitting phase boundaries and adsorption datajs, 16]. However, while 
this model is closest to the one we analyze, the two differ in a number of points. For example, 
in contrast to our model the PP, HH and PH interactions within the Karlstrom model are 
not identical. With this caveat in mind, this set of fi° H , fi P , fi SP and ix° SH can be used to 
illustrate the behavior of e e // for a "PEO-like" case. In the Karlstrom model the standard 
chemical potential of the two-states in the bulk is given by fi® = Ui — RT In ft where Uj is 
the internal energy and ft is a degeneracy factor. The values per mole in KJ are fi° P = 
and /i p = 5.086 — RT In 8 [8| . For the adsorbed species /ig; — ^ — Ae" d where the adsorption 
energy per mole Aef at methylated silica are Aef = 0.825ATJ and Ae^f = 1.Q25KJ [16f. 
A rough idea concerning the physical consequences of the Karlstrom model may be gained 
upon comparing the e e // of a hypothetical P homopolymer, corresponding to PEO modeled 
as hydrophilic one-state polymer, to the annealed PH homopolymer model of PEO. Such 
comparison shows that e e ff of the two-state chain is shifted upwards by roughly a factor of 
two. 



IV. SIMULATION MODEL 



We simulate a two-state polymer terminally anchored to a planar surface. Following 



Baumgartner 



321 ] the polymer is modeled as freely jointed chain comprising N Lennard- 



Jones (LJ) particles. The monomers within this bead-spring model interact via a LJ poten- 
tial 



V, 



L.I 



12 



(25) 



where r. 



'.i 



|r; — r j | is the distance between monomers i and j such that \i — j\ > 2 



and Ti is the position vector of the zth monomer, e specifies the depth of the potential 
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minimum at r = 2 1 / 6 a and a, the collision diameter, is the separation for which Vlj = 0. 
This potential exhibits a soft-core steric repulsion at r < a, and steeply decaying attraction 
for r > a. The monomers exist in two interconverting states, P and H. It is thus necessary 
to specify three LJ potentials, corresponding to the interactions between PP, HH, and PH 
monomer pairs, involving altogether six parameters. In the following all three potentials are 
characterized by the same a and we will thus express all distances in these units. The three 
remaining parameters, ? = eppkT; epn&T; e^nkT determine the strength of interactions 
between the PP, PH, and HH monomers. For the simulation of self avoiding chains we 
set epH = ?hh — £pp — 0.20 so that the monomer- monomer interactions are dominated by 
the excluded volume and contain effectively no attractive contribution. For simulations of 
ideal polymers <Tph = ?hh = ?pp = 0. In addition, the bulk HP states are characterized 
by standard chemical potentials fi° H = kTAe and jjp — with Ae > 0. Accordingly 
A// = AekT specifies the difference in standard chemical potential A// = fi° H — fi° P between 
noninteracting H and P monomers. All monomer pairs, except nearest-neighbors, interact 
via LJ potentials. Nearest-neighbor monomers along the chain are constrained to 0.7 < 
(r i+ i — Ti) 2 < 2.0. Separations outside this range incur an infinite energy penalty thus 
ensuring connectivity. 

We have used two monomer-surface interaction potentials, t/ wa ii, both specified in terms 
of z, the distance between the monomer and the surface. One is the contact potential 



oo z < 



wall 



-kTAef < z < 1 (26) 
z > 1 

where Ae°; d > is the depth of the square well adjacent to the surface as experienced 
by species i = P,H. We will mostly focus on the contact potential because it allows for 
unambiguous definition of adsorption i.e., a monomer having its center within the slab 
< z < 1.0 is adsorbed. The second, more physical, "10-4" wall potential 



C/ wall = kTAe" 



ad 



(27) 



as obtained by integrating the LJ potentials between the top monolayer atoms of the sub- 
strate and a monomer at z. Here again Ae" d > is a tuning parameter characterizing the 
strength of the attraction [3]. When using the 10-4 potential the criterion for adsorption 
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is somewhat arbitrary and we define monomers with z < 2 as adsorbed [3]. Note however 
that monomers within the adsorption slab experience a z dependent U wa x\. In both cases the 
values of Aef d for the P and H states, Ae^f or Ae P d , differ favoring the surface H state i.e., 
Ae p d < Ae^f. For the contact surface potentials we thus have 

fi° SH /kT = fi° H - Ae# = Ae - Ae# and fi° SP /kT = fi° P - Ae a P d = -Ae a P d . (28) 

For the 10-4 potentials fi° SH and dependent. However, due to the fast decay of the 

potential /j° sh and fi° SP as given by ( |28i) provide, as we shall see, a reasonable approximation. 

The simulations involves chains comprising N = 64, 128, 256 monomers. At each Monte 
Carlo step (MCs) we shift the position of every monomer in the chain and update its HP 
state using the Metropolis algorithm. This procedure is thus repeated N times per MCs. 
For each set of parameters the simulation involves 2 x 10 7 MCs. The system is equilibrated 
during million MCs. The remaining 1.9 x 10 7 MCs are grouped into sets of 10 4 MCs whose 
configurational and PH characteristics are averaged for analysis. In a conformational up- 
date of the monomer position, a monomer is chosen randomly, its position is shifted by 
a sufficiently small distance and the resulting energy difference kTAE, accounting for LJ 
interactions, is calculated. When AE < 0, the operation is accepted and we proceed to the 
next monomer movement. On the other hand when AE > 0, the operation is accepted with 
the probability exp(— AE). The PH interconversion updates are implemented in two stages 
corresponding to the bulk equilibration and its modification by the surface potential. In the 
"bulk stage" a monomer is randomly chosen and its state is updated with the probability 

p(P -> H) = exp(-Ae) 

(29) 

p(H -> P) =1 

where p(P — > H) and p(H — > P) denote respectively the transition probabilities from P to H, 
and from H to P. The detailed balance condition for the bulk equilibrium is p(P — > H.)N PP = 
p(ll^P)N BH yields 

^£ = ^ = exp(-Ae)^. (30) 

This procedure is the counterpart of the trial motion in the conformational steps. The 
PH states are subsequently updated, in z dependent fashion, allowing for the effect of the 
surface potential on the PH equilibrium. For the contact potential, monomers with z = 1 
are converted following the transition probabilities 
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P{SP -> SH) = exp(Ae^ - Ae^) 
P(SH -> SP) =1 

where P(SH — > SP) = 1 because Ae a p < Ae# . The corresponding detailed balance condition 
for equilibrium at the surface, p(P -> H)P(SP -> SH)N SP = p(H -> P)P(SH -> SP)N SH 
yields 

A^£ = = exp(A ^ _ Ae _ Ae a d) = ^ (32) 

In contrast, monomers at z > 1 experience no surface effect and is accordingly unmod- 
ified. Similar procedure is used for the 10-4 wall potential with the distinction that the z 
dependent surface transition probabilities occur, in principle at all z values with subsequent 
effect on Kb as well as on K$- However, since the 10-4 potential decays fast the effect on 
the bulk equilibrium is negligible. 



V. SIMULATION RESULTS 



As opposed to the case of quenched copolymers, whose composition is fixed, the H 
fraction of our model two-state polymers depends on three tuning parameters Ae, Aepf 
and Ae a P d . In the following Ae a P d is fixed while Ae and Ae ^ are varied. Our results are 
mostly concerned with chains exhibiting self avoiding random walk (SAW) statistics and 
interacting with the surface via contact potentials. In addition we will comment briefly 
on the behavior of ideal chains exhibiting random walk (RW) statistics and the effect of 
10-4 monomer-surface potentials. Importantly, variation of Ae%f affects Nsh (Figure 2), 
N H (Figure 3) and A^" 5 (Figure 4). As seen from Fig 2, N SH /N S = K s /(1 + K s ) and 
N BH /N B = K B /(1 + Kb) in agreement with equations (jlj) and (j3J). To proceed further it 
is first necessary to determine e^y. As noted earlier, curves of Ng/N vs. e e ff for different 
iV values intersect at e e ff = e c e ^ provided that $ = 1/2. Alternatively, curves of N$/N 
vs. Aejf will intersect at Ae^ for families of curves of identical Ae, Acp*, but different iV 
values. Within the accuracy of our data the curves do intersect at a single point (Figure 5) 
thus lending support to the assumption that the cross-over exponent $ = 1/2 describes the 
adsorption behavior of two-state polymers. Furthermore, as obtained by substituting 
Ae^ and the corresponding Ae, Ae a p values into (jHJ) is a constant (Figure 6) as expected 
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within the picture of "homopolymer-like" adsorption. Having obtained e c e ^ we find that the 
Ns vs Ae ^ data collapses onto a universal curve (Figure 7) upon plotting Ns/N 1 / 2 vs tN 1 / 2 
where 

T _ kT ln ex P ( Ae ^ ~ Ae ) + ex P ( Ae p) x f33 A 
11 1 + exp (-Ae) ' 1 j 

Furthermore, the N B /N data agrees with equation (12B as adapted to the simulation model 
(Figure 3) 

Nh „ exp(Ae + Ae a P d -Ae a H d ) _ exp(Ae) 

AT ~ exp(Ae + Ae<^ - Ae^) + 1 1 J exp(Ae) + 1 ' 1 ' 

The above results concern SAW and contact potentials. The e c e ^ values for RW are lower 

than those of SAW chains and depend on the wall potential (Figure 6). However, since 

$ = 1/2 applies to SAW as well as RW chains, the Ns/N 1 / 2 vs tN 1 / 2 scaling with the 

proper choice of r is obtained in both cases for adsorption due to contact potentials (Figure 

7-8). 

The adsorption behavior of homopolymers experiencing 10-4 potential and contact po- 
tential is indistinguishable [3]. Our discussion suggests that the adsorption of an annealed 
copolymer is analogous to the adsorption of a homopolymer with e e jf given by (jHJ). Ac- 
cordingly we expect the adsorption of annealed homopolymers to be insensitive to the 
choice of the potential and of the z values used to define adsorbed monomers. Within 
the accuracy of our simulation, this is indeed the case for the adsorption of RW annealed 
copolymers where the data is collapsed by plotting Ns/N 1 ^ 2 vs tN 1 ^ 2 (Figure 9). It is how- 
ever important to note that for this case the PH equilibrium is z dependent. Accordingly 
Nsh/Ns = K s /(1 + K s ) and N BH /N B = K B /{1 + K B ) as specified by flU and ([32} only 
capture the qualitative features of the PH equilibrium but do not yield quantitative agree- 
ment because fj^ H and n° SP as approximated via (f2"g|) do not allow for the z dependence 
of t/waii- A quantitative agreement is achieved upon dividing Ae# or Ae P d by 1.5 allowing 
presumably for the average U wa \\ experienced within the z < 2 region (Figure 10). 

VI. DISCUSSION 

The adsorption behavior of annealed copolymers may conceivably differ with the details 
of the model. Our discussion concerned the particular case of non-cooperative two-state 
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polymers involving unimolecular PH interconversion. For this case, our results suggest that 
the annealed copolymers adsorb as homopolymers once the appropriate effective adsorption 
energy per monomer at the surface, e e //, is identified. This e e ff accounts for the average 
adsorption energy as well as the mixing entropy of the surface PH monomers at equilib- 
rium. The adsorption of two-state polymers is described by the cross-over exponent <3> of 
homopolymer adsorption. In this respect, their behavior is identical to that of quenched 
copolymers. However, as opposed to quenched copolymers, where e e // as obtained from 
the MW annealed approximation is determined by the fixed average H/ P fraction 22|, |25| , 
e e ff in the annealed case depends on the PH equilibrium constants in the bulk and at the 
surface, Kb and K$- In other words e e ff varies with the H/P ratio which, in contrast to the 
quenched case, depends on the temperature, the adsorption energies and the bulk standard 
chemical potentials. This picture allows to calculate the number of adsorbed monomers Ns, 
the H/P ratio in the bulk and at the surface as well as the overall H/P ratio. The calculated 
values are in good agreement with the simulation results concerning SAW and RW experi- 
encing contact potentials. Simulation of homopolymer chains suggest that the form of the 
wall- monomer potential does not affect $ and the scaling behavior of N$- Since annealed 
copolymers adsorb as homopolymers characterized by e e ff one expects similar insensitivity 
to the details of the potential. In contrast, the H/P ratio is sensitive to the range of the 
wall potential. Both features are apparent from our results concerning the adsorption of RW 
annealed copolymers subject to a 10-4 wall potential. As noted earlier, our analysis utilizes a 
simplified version of the Karlstrom model. In contrast to this model all monomer-monomer 
interactions are identical. New features may well occur upon introducing attraction between 
hydrophobic H monomers. In the case of free chain in the bulk such interaction can lead to 



phase separation 18] and to polymer collapse 



33l |. The preferential adsorption of H monomers 



together with HH attraction may lead to adsorption induced collapse. 

Our discussion thus far focused on the comparison between simulation results of annealed 
and quenched copolymers. The comparison of the experimental situation brings up a second 
issue. The experimental realization of quenched copolymers is clear and studying the ad- 
sorption of the corresponding homopolymers allows to deduce the adsorption energies of the 
different monomers involved. The situation is more difficult in the case of annealed copoly- 
mer as encountered in the modeling of NWSP. The corresponding one state homopolymers 
do not exist. Furthermore, direct experimental evidence concerning the molecular identity 
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of the interconverting states is yet to emerge. As a result the interaction parameters are 
obtained indirectly by model dependent fitting of experimental data. With this in mind 
it is important to recall the experimental evidence for two-state models. Two observations 
are of special interest. One is their ability to predict the qualitative features of the phase 

diagram and fit experimentally observed phase boundaries [7|, y,y,Q,Q,Q,Q,Q. The 

second concerns the stretching of PEO chains in atomic force microscopy experiments. This 
reveals that different force laws characterize the strong stretching regime in water and in 
hexadecane. In particular, the chain extension in water exhibits a plateau characteristic of 
two-state polymers|34j]. 



Acknowledgement 1 N.Y. benefited from fellowship No. 7662 from the Japan Society for 
the Promotion of Science (JSPS) as well as from the hospitality of the theory group of the 
ILL. Part of the numerical calculations in this work was carried out on Altix3700 BX2 at 
YITP in Kyoto University. 



[1] Napper, D. H. Polymeric Stabilization of Colloidal Dispersions; Academic Press: London, 
1983. 

[2] Eisenriegler, E. Polymers Near Surfaces: Conformation Properties and Relation to Critical 

Phenomena; World Scientific Publishing Company 1993 
[3] Molyneux, P. Ed. Water-Soluble Synthetic Polymers: Properties and Behavior; Volume I &II 

CRC Press 1984 

[4] Glass J.E. Ed, Water-Soluble Polymers: Beauty with Performance; Oxford University Press, 
USA 1986 

[5] Glass, J.E. Ed., Hydrophilic Polymers: Performance with Environmental Acceptance; Oxford 
University Press 1996 

[6] Williams, P. Ed. Handbook of Industrial Water Soluble Polymers; Wiley-Blackwell NY, 2007 
[7] Goldstein, R. J. Chem. Phys. 1984, 80, 5340-5341. 
[8] Karlstrom, G. J. Phys. Chem. 1985, 89, 4962-4964. 
[9] Matsuyama, A.; Tanaka, F. Phys. Rev. Lett. 1990, 65, 341-344. 
[10] Veytsman, B. J. Phys. Chem. 1990, 94, 8499-8500. 



16 



[11 

[12 
[13 
[14 
[15 
[16 

fir 

[18 

[19 
[20 
[21 
[22 
[23 
[24 
[25 

[26 

[27; 

[28 
[29 



de Gennes, P.G. CR I'Academie. Sci., Ser. II 1991 313, 1117-1122. 

Bekiranov, S.; Bruinsma, R.; Pincus, P. Europhys. Lett. 1993, 24, 183-188. 

Jeppesen, C; Kremer, K. Eurphys. Lett. 1996, 34, 563-568. 

Bekiranov, S.; Bruinsma, R.; Pincus, P. Phys. Rev. E 1997, 55, 577-585. 

Dormidontova, E. Macromolecules 2002, 35, 987-1001. 

Linse P; Bjorling M. Macromolecules 1991, 24, 6700-6711 

Eisenriegler, E.; Kremer, K.; Binder, K.; J. Chem. Phys. 1982, 77, 6296 

Metzger, S.; Miiller, M.; Binder, K.; Baschnagel, J. Macromol. Theory Simul. 2002, 11, 985- 

995 

Luo MB J. Chem. Phys. 2008, 128, 044912 

Bouchaud, E; Daoud, M J. Phys. (Orsay) 1987, 37, 1991-2000, 

Sumithra K; Baumgaertner, A. J. Chem. Phys 1999, 10, 2727-2731 

Moghaddam, MS; Whittington SG J. Phys. A: Math. Gen, 2002, 35, 33-42 

Soteros, CE ; Whittington SG J. Phys. A: Math. Gen, 2004, 37, R279-R325 

Ziebarth, J D.; Wang, Y; Polotsky,A; Luo M. Macromolecules 2007, 40, 3498-3504 

Bhattacharya, S.; Hsu,H.-P.; Milchev, A.;Rostiashvili,V. G.; Vilgis, T. A. Macromolecules 

2008, 41, 2920-2930 

Grosberg, A. Yu.; Khokhlov, A. R. Statistical Physics of Macromolecules; American Institute 
of Physics, New York, 1994. 

de Gennes, P.-G. Scaling Concepts in Polymer Physics; Cornell, University Press: Ithaca, NY, 
1979. 

Rubinstein, M.; Colby, R. H. Polymer Physics; Oxford University Press: Oxford, U.K., 2003. 
A fixed overall H fraction is not gauranteed in the MW annealed model especially for large e 



23(]. 



[30] Our discussion focused on the MW annealed approximation because of its simplicity and 



proven utility in the analysis of simulation data 



24 



25]. We should add that the more general 



Morita approximation 23J is also obtainable from equation (|10p upon choosing rfp/kT 
H° SP /kT = Xp — ln(l — p), —fi° H /kT = A(l — p) +lnp and —fjP SH /kT = €mw — hip, n Q SH /kT 
£mw +A(1 -p) + lnp. 

[31] Hsu, H.-P, Nadler, W.; Grassberger P J. Phys. A: Math. Gen, 2005, 38, 775-806 

[32] Baumgartner, A. J. Chem. Phys. 1980, 72, 871-879. 

17 



[33] Yoshinaga, N. ; Bicout, D.J.; Kats E.I.; Halperin, A. Macromolecules 2007, 40, 2201 
[34] Oesterhelt, F., Rief M., Gaub, H., New J. Phys. 1999, 1, 6. 



18 



(A) 



FIG. 1: Two configuration of an adsorbed two-state polymer undergoing P(full circles) H(empty 
circles) interconversion in the bulk and at the surface. 
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FIG. 2: Nbh/Nb and Nsh/Ns as a function of a contact potential Ae— Aejf for various fj^—iufp = 
kTAe and the corresponding curves given by ([30]) and (|3"2"]h 
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FIG. 4: Ns /N vs Ae^f of a contact potential for self avoiding chains with various [jP h — /ip = kTAe 
and N. 
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FIG. 6: Ae c H and the corresponding e^,, (as obtained by substituting Ae c H into (jHJ)) vs Ae. 
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FIG. 7: Ns/N 1 / 2 vs riV 1 / 2 for SAW experiencing a contact potential. 
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FIG. 8: Ns/N 1 / 2 vs riV 1 / 2 for RW experiencing a contact potential. 



22 





FIG. 10: N BH /N B and N SH /N S as a function of a 10-4 potential Ae - Aeff/1.5 for various 
fj,jj — (jJp = kTAe and the corresponding curves given by ([30]) and ([32]) , 
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